Meissner state in finite superconducting cylinders with uniform applied magnetic field 



OS 
Os 
On 



O 



o 
o 

c3 



i 

-a 
c 
o 
o 



> 

in 

O 
O 

On 

o> 



-a 

c 
o 
o 



s 



F. M. Araujo-Moreira 

Grupo de Supercondutividade e Magnetismo, Departamento de Fisica, 

Universidade Federal de Sao Carlos, Caixa Postal 616, Sao Carlos SP, 13565-905 Brazil 

C. Navau and A. Sanchez 
Grup d'Electromagnetisme, Departament de Fisica, Universitat Autdnoma Barcelona 08193 Bellaterra, Barcelona, Catalonia, 

Spain 

We study the magnetic response of superconductors in the presence of low values of a uniform ap- 
plied magnetic field. We report measurements of DC magnetization and AC magnetic susceptibility 
performed on niobium cylinders of different length-to-radius ratios, which show a dramatic enhance 
of the initial magnetization for thin samples, due to the demagnetizing effects. The experimental 
results are analyzed by applying a model that calculates the magnetic response of the superconduc- 
tor, taking into account the effects of the demagnetizing fields. We use the results of magnetization 
and current and field distributions of perfectly diamagnetic cylinders to discuss the physics of the 
demagnetizing effects in the Meissner state of type-II superconductors. 



I. INTRODUCTION 

Important intrinsic parameters of superconductors, 
such as the lower critical field H c ± and the critical cur- 
rent density J c , are experimentally obtained by measur- 
ing their response to an applied magnetic field. The pro- 
cedures to obtain these parameters often rely on theo- 
retical approaches developed for infinite samples. When 
considering realistic finite-size superconductors, impor- 
tant complications arise, so that these methods fail. Even 
in the case of a uniform field applied to a finite super- 
conductor, in general it is not easy to extract informa- 
tion about the intrinsic parameters of the superconduc- 
tor since the magnetic response may be strongly depen- 
dent on its shape. Demagnetizing effects appearing in 
finite samples make the internal magnetic field H— B//xo 
in the sample different than the applied one, H a . The ex- 
act relation between both magnetic vectors is in general 
unknown. As a consequence of this indetermimation of 
the local magnetic field in the sample, the estimations of 
H c i and other parameters are very complicated. Some- 
times, this problem is treated by considering a constant 
demagnetizing factor N, so that the magnetic field in the 
sample volume H is assumed to be related to the applied 
magnetic field H a by: 



H = H n 



NM. 



(1) 



This procedure is correct only for ellipsoidal-shape sam- 
ples and when the external magnetic field is parallel to 
one of its principal axis JD . In all other cases this equa- 
tion is not valid and it becomes very hard to find a sim- 
ple relation between the magnetic vectors H and H a . 
Moreover, in general the field H is related with H Q in a 
different way from point to point in the superconductor. 
Recently, there have been important theoretical ad- 
vances in treating the magnetic response of finite super- 
conducting samples. The magnetic response of finite su- 



perconductors in the critical state including demagne- 
tizing effects have been recently calculated for strips 0] 
and cylinders |3|-|5[] , following previous works on very thin 
strips and disks JTJ]. These works deal basically with 
current and field distributions in superconductors in the 
mixed state, with critical-state supercurrents penetrat- 
ing into the bulk of the material. Besides, Chen et al. Q 
have calculated demagnetizing factors for cylinders as a 
function of the length to diameter ratio and for different 
values of the susceptibility x (including \ — — 1)- 

However, to our knowledge, there has not been done 
a systematic study involving the comparison between 
experiments and theoretical data of the Meissner state 
in superconducting cylinders. This is equivalent to ask 
which are the currents that completely shield a cylindri- 
cal volume and the resulting magnetization. 

In this work we systematically investigate the magnetic 
response of superconducting cylinders in the complete 
shielding state, by quantitatively studying the effect of 
demagnetizing fields in their Meissner response. This pa- 
per is organized as follows. In Section II we discuss the 
experimental setup and measured samples. In Section 
III we present the experimental results obtained from 
DC magnetization and AC magnetic susceptibility tech- 
niques, performed on niobium cylindrical samples. We 
introduce in Section IV a new approach to study the mag- 
netic response of completely diamagnetic finite cylinders. 
This model allows the interpretation and understanding 
of our experimental results, which is discussed in Section 
V. The model enables us to calculate the surface current 
distributions resulting from the magnetic shielding of the 
cylinder, as well as the magnetic fields created in the ex- 
terior of the samples by the induced supercurrents. The 
results are discussed in Section VI. Finally, we summarize 
our conclusions in Section VII. 



II. SAMPLES AND EXPERIMENTAL SETUP 

We have performed the magnetic characterization of 
ten cylindrical niobium polycristalline samples with dif- 
ferent values for the ratio L/R, where L and R are the 
length and the radius of the sample, respectively (see Ta- 
ble I). They were obtained from two different pieces of 
brut niobium. After machined, they were cleaned in ul- 
trasound bath, and later by using a HCI-HNO3 solution. 

We have studied two families (with five samples in each 
one) of cylinders with nominal diameters of 1.94 mm and 
2.86 mm, respectively. Samples with the same nominal 
value for the diameter were cut from the same piece, 
in cylinders of different lengths. The variation between 
nominal and real values of the diameter is smaller than 
4%. 

Before performing magnetic measurements, we have 
determined the quality of all samples through x-ray 
diffraction (XRD) and scanning cletron microscopy 
(SEM) analysis. For the XRD we use a SIEMENS 
D/5000 difractometer, and for the SEM, we use a JEOL 
JSM-5800/LV microscope. 

The magnetic characterization was performed through 
both the magnetization as a function of the external DC 
magnetic field, M(H a ), and the complex AC magnetic 
susceptibility, as a function of the absolute temperature, 
x(T). To perform those experiments we have used a 
Quantum Design MPMS5 SQUID magnetometer able to 
operate in the ranges 2 K < T < 400 K, 0.1 A/m < 
h < 300 A/m, and 1 < / < 1000 Hz, where T, h and / 
are the absolute temperature, and the amplitude and the 
frequency of the AC magnetic field, respectively. In all 
cases, to avoid trapped magnetic flux, samples were zero- 
field-cooled (ZFC) before each experiment. The magnetic 
field was always applied along the axis of the cylinders. 



III. EXPERIMENTAL RESULTS 

The XRD and SEM analysis confirmed the high qual- 
ity of the niobium samples, as verified by the absence 
of impurities and the low density of grain boundaries. 
Measurements of the complex AC magnetic susceptibil- 
ity, x(T)=x'(T)+ x"(T), were performed with the pa- 
rameters h = 80 A/m, and / = 100 Hz, and are shown in 
Fig. 1. These experiments also confirm the quality of our 
polycristals, by showing a critical temperature Tq = 9.2 
K, and sharp superconductiong transitions, with typical 
width of 0.2 K. This value is considered excellent for a 
polycristalline sample. These measurements also point 
out the strong shape effect on x{T)- As can be seen in 
Fig. 1, the lower the ratio L/R is, the larger the shape 
effect is, evidenced for larger values of the modulus of 
X'(T). ' 

As mentioned above, we have also verified the mag- 
netic behavior of the cylinders by measuring isothermal 
M(H a ) curves. The measurements were performed at 



T = 8K. We show these results in Fig. 2. There, we can 
see again the strong shape effect on the obtained curves. 
The lower the ratio L/R, the higher the value of the ini- 
tial slope M/H . Also, we can see the shape effect on 
the point at which the magnetization curve departs from 
a straight line, which corresponds to the position where 
the flux starts to penetrate in the superconductor. Thus, 
flux penetration starts at lower fields for large values of 
L/R, as expected. 

As it was mentioned in the last section, samples with 
the same nominal diameter were cut from the same piece, 
in cylinders with different lengths. Since the variation be- 
tween nominal and real values of the diameter is smaller 
than 4%, only shape effects will be responsible for the 
observed difference in their magnetic response. 



IV. DESCRIPTION OF THE MODEL 

In order to study the Meissner state and analyze the 
experimental data, we have developed a model to calcu- 
late the surface currents that shield any axially symmet- 
ric applied magnetic field inside a cylindrical material. 
This model could be applied, in principle, to both type-I 
and type-II superconductors below H c and H c \, respec- 
tively, and also to good conductors in a high frequency 
applied magnetic field. 

It is well known that supercurrents are induced in the 
superconductors in response to variations of the applied 
magnetic field. For zero-field cooled type-II superconduc- 
tors in the Meissner state, in order to minimize the mag- 
netic energy, these supercurrents completely shield the 
applied magnetic field inside the superconductor. This 
shielding occurs over the whole sample volume except for 
a thin surface shell of thickness A (the London penetra- 
tion depth), where supercurrents flow. Our model simu- 
lates this process. It will enable us to calculate the cur- 
rent distribution that shields the applied magnetic field, 
by finding at the currents that minimize the magnetic 
energy in the system. For simplicity, we will assume that 
A is negligible. In this approximation, the calculated su- 
percurrents flow only in the cylinder surfaces. 

We consider a cylindrical type-II superconductor of 
radius R and length L with its axis in the z direc- 
tion, in the presence of a uniform applied magnetic field 
H a = H a z. We use cylindrical coordinates (p, ip, z). Ow- 
ing to the symmetry of the problem all shielding currents 
will have azimuthal direction. We divide the supercon- 
ducting cylinder surfaces into a series of concentric circu- 
lar circuits in which currents can flow, with no limitation 
in the value of current. We consider n of such circuits 
in each one of the cylinders ends and m circuits in the 
lateral surface (see Fig. 3). These linear circuits are 
indexed, as seen in Fig. 3, from i = I to i = 2n + m — 1. 

The method to obtain the current profile is the follow- 
ing. We start with a zero-field cooled (ZFC) supercon- 
ductor and set an applied field H a . The magnetic flux 



that threads the area closed by the i— circuit due to the 
external field is: 



*? = HonplHa, 



(2) 



Pi being the radius of the i— circuit. The magnetic flux 
contributes with an energy that has to be counteracted 
by induced surface currents. A step of current AI that 
is set in some j— circuit requires an energy: 



E 



. - \l 3 ^i? 



(3) 



being Lj the self-inductance of the j— circuit, while it 
contributes to reduce the energy by a factor (A7)$". So, 
after an external field is applied, we seek for the circuit 
that decreases the energy the most and set a current step 
AI there. The same criterion is used next to choose either 
to increase the current at the same circuit by another 
step AI, or instead put AI at some new circuit. In the 
last case (and in the general case when they are curents 
circulating in many circuits of the sample), the energy 
cost to set a current step has an extra term coming from 
the mutual inductances of all other currents: 






(4) 



The mutual inductances M^ between the i and j lin- 
ear circuits are calculated using the Neumann formulas 
(see, for example, Ref. [8]). To avoid diverging self- 
inductances, we have used a cut-off and calculated the 
self-inductance of a circuit with radius p from the mu- 
tual inductance between the considered circuit and one 
with the same current at a radial position p + e. An ap- 
propriate choice for e has been found to be e = 0.78Ai?, 
where AR is the radial separation between two consecu- 
tive circuits in the end faces (AR — R/n). 

The minimun energy correspoding to a given value of 
the applied magnetic field will be reached when it be- 
comes impossible to further decrease the magnetic en- 
ergy by setting extra current steps. From the existing 
current profile, we can easily obtain the different mag- 
nitudes we are interested in. The magnetization, which 
has only axial direction M z , is calculated by using: 



M 7 . 



1 



R 2 L 



E 1 ^ 2 



(5) 



The magnetic induction B could be computed from the 
current profiles using the Biot-Savart formulas. However, 
we use a simpler way, which allows the calculation of B 
from the flux. In our model, the total magnetic flux that 
threads any circular circuit (not necessarily those in the 
surfaces of the superconductor) can be easily calculated 
as: 



$(p,z) = $ a (p,z) + &(p,z), 



(6) 



where the internal flux that threads a j— circuit due to 
all the currents is <f>* = J^k Mjklk (the term of the self- 
inductances is also included) . Then, the axial component 
of the total magnetic induction is simply: 



B z (p,z) = 



<P(p + AR,z)-<S>{p,z) 
tt((p + AR)2-p 2 ) ' 



(7) 



The radial component B r (p, z) is calculated from the val- 
ues of the axial one, imposing the condition that the di- 
vergence of B equals zero. 

Finally, if the external magnetic field is further in- 
creased, the same procedure starts again from the present 
distribution of current. 

The values of n and m have been chosen sufficiently 
large so that the results are independent of their partic- 
ular value. Typical values are n ■ m ~ 7500 — 10000. The 
computation time for a initial curve in a personal Digital 
Workstation takes few minutes for any L/R ratio. 

We would like to remark that, with the method de- 
scribed here, we obtain the current distribution, field 
profiles, magnetization, and all the other results without 
using any free parameter. Only the direction of currents 
has to be known, which is straightfordward in this geom- 
etry. 



V. COMPARISON OF EXPERIMENTAL AND 
THEORETICAL DATA 

In Fig. 4 we compare the experimental values of the 
initial magnetic susceptibility ^; n ; obtained from both 
DC and AC measurements, with those calculated from 
our model, for different values of L/R. The calculated 
values of the initial slope are only function of the length- 
to-radius ratio L/R. The agreement with the experi- 
mental data is then very satisfactory, confirming the va- 
lidity of our theoretical approach. Both experimental 
and calculated data indicate a strong increase in the ab- 
solute value of the initial susceptibility when decreasing 
the cylinder aspect ratio. When the sample is very large 
(for L > 10-R), the initial susceptibility X'mi approaches 
the value predicted for infinite samples, ^ini = — 1- For 
shorter samples, the magnetization gets larger in magni- 
tude for the same value of the applied magnetic field. We 
find that the general behavior of the dependence of xi n i 
on L/R is well described (with a departure of less than 
1.5 % from our calculations) by the approximate formula 
given by Brandt (pf: 



-ttR 2 L 



s 



R s 



AR C 



tanh 



1.27 — In 
2R 



1 



2R 
T 



(8) 



The values of the initial slope calculated from our model 
are also compatible with those given by Chen et al. Q 
for cylinders with x = — 1 with a maximum deviation of 
1.0%. 



VI. DISCUSSION 

A. Effect of the demagnetizing fields 

The experimental and theoretical results can be un- 
derstood from analyzing the effect of the demagnetizing 
fields. In an irrealistic infinitely long sample, if the field 
is applied along its axis, shielding currents will flow only 
in the lateral surface, with a constant value along the 
cylinder. This creates a spatially uniform magnetic field 
over the superconductor, which makes B=0 inside. In 
finite samples, however, in the top and bottom end sur- 
faces the tangencial magnetic field is not continuous and 
shielding currents are also there induced. Hence, these 
currents create an extra non-uniform magnetic field over 
the lateral surface of the cylinder, so the currents flowing 
in this surface will not have a constant value. It is easily 
seen (by a simple examination of the magnetic field cre- 
ated by each single current loop) that the effect that the 
demagnetizing fields produce in the lateral surface region 
is to enhance the local magnetic field. In this case, its 
total value is larger than the actual value of H a . As a 
result, higher values of current are necessary to shield the 
applied magnetic field, yielding to larger values of both 
the magnetization and the magnetic susceptibility. It is 
clear that, the thinner the sample is, the larger this effect 
occurs, which explains the increase of \x\ when decreasing 
L/R. 



B. Field and current profiles 

The previous discussions can be illustrated by studying 
the distribution of the magnetic field in finite cylinders. 
In Fig. 5 we show the calculated total magnetic induc- 
tion B for cylinders with three different values of L/R. 
The displayed lines indicate the direction of B (tangen- 
tial to the lines at each point), although their densities 
do not reflect in general the field strength. These results 
show that only in the case of the largest L/R ratio the 
magnetic field in the lateral surface can be considered 
as basically having a constant (axial) direction over the 
cylinder length. In the other two cases, the magnetic field 
loses its main axial direction, gradually bending towards 
the axis at both cylinder ends. 

In Fig. 6 we plot the calculated surface current den- 
sity corresponding to the L/R = 1 case. In both top and 
bottom ends of the cylinder, the strength of the shielding 
supercurrents flowing in the azimuthal direction gradu- 
ally grows (in absolute value) from a zero value over the 
axis to a diverging behavior when reaching the cylinder 
edge (this divergence is smoothed because of our dis- 
cretization; in actual experiments the nonzero value of 
A makes also smooth values) . The currents in the lateral 
surface are also stronger at the edges and their intensity 
decrease towards the center of the sample, where they 
have a roughly constant value over the plateau shown in 



Fig. 6. The extension of this plateau (if defined as the 
region where the surface current differs less than 5% with 
respect to the minimum value, located at the center) is 
of about 70 % of the total cylinder length, for L/R = 10. 
This percentage decreases to 46% and 36% for L/R = 1 
and L/R = 0.2, respectively (for clarity, Fig. 6 does not 
show the data for the cases L/R = 0.2 and 10). This is in 
correspondence with the regions, shown in Fig. 5, where 
the magnetic field is almost constant. Besides, our calcu- 
lations show that the relative contribution of the currents 
in the end surfaces to the magnetization increases by de- 
creasing the cylinder thickness. The contributions of the 
lateral and the end surfaces to the total magnetization 
are depicted in Fig. 7. These results show that, whereas 
in a long sample (L/R — 10) about 94% of the contri- 
bution to the superconductor magnetization comes from 
the lateral surface, this percentage decreases to about 
72% and 43%, for L/R = 1 and L/R = 0.2, respectively. 



C. Remarks about the generality of the results 

In the M(H a ) experimental data, Fig. 2, the curves 
for different L/R show a systematic behavior in the first 
portion of the initial curve (region of interest). There, 
the shielding should be perfect. Nevertheless, the trend 
is not so clear for larger values of the applied magnetic 
field. In the middle part of the loop, it is known that bulk 
supercurrents penetrate into the superconductor, which 
goes from the Meissner state to the mixed state. Within 
the critical-state model framework, the magnetization in 
the mixed state is known to be dependent on the partic- 
ular J C (B) dependence of the samples H. This depends 
itself on factors such as the detailed microstructure of 
the sample. Thus, samples with different microstructures 
may have a difference distribution of pinning strenghts 
and locations, which influences the critical current and 
the magnetization. The samples we have measured were 
cut from two different cylinders, which explains why the 
M(H a ) loops for intermediate values of H a are different. 
The systematic behavior observed in the initial suscepti- 
bility despite using samples from different original pieces, 
supports the generality of our results for any diamagnetic 
cylinder. 



VII. CONCLUSIONS 

In this work we propose a new model based on energy 
minimization which allows the calculation of the mag- 
netic response for perfectly diamagnetic cylinders of any 
size with high precision. The only assumption we have 
considered is that the magnetic penetration depth, A, is 
negligible. We have experimentally verified the validity 
of this model by measuring the low-field magnetic re- 
sponse of niobium cylinders with different values for the 



lcngth-to-radius ratio. We have demonstrated both ex- 
perimentally and theoretically that the value of the ini- 
tial susceptibility of zero-field cooled type-II supercon- 
ductor cylinders is a function of only the sample aspect 
ratio, and calculated its value for a wide range of sam- 
ple dimensions. The model is sufficiently general to be 
adapted to other geometries and also to non-uniform ax- 
ially symmetric applied fields. These results may help 
to discriminate whether some of the effects recently ob- 
served in the study of thin film high-T c superconductors 
are due to intrinsic causes or instead have an extrinsic 
origin associated with sample size effects. 
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FIG. 1. Real part of the AC magnetic susceptibility as a 
function of temperature, x'(T), performed with the parame- 
ters h = 80A/m, and / = 100Hz. Samples identifications are 
listed in Table 1. 
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TABLE I. Values of length (L) and radius (R) for the ten 
Nb samples studied 
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FIG. 2. DC Magnetization loops as a function of the ap- 
plied magnetic field, M(H a ), at T — 8 K. Samples identifica- 
tions are listed in Table 1. 



FIG. 3. Sketch of the discretization of the cylinder used in 
the model. 



FIG. 4. Initial susceptibility for diamagnetic cylinders of 
radius R and length L as a function of L/R. Squares and 
circles are experimental data from AC and DC experiments 
performed on niobium cylinders, respectively. The line cor- 
responds to our model calculation. Vertical error bars are 
smaller than symbol size. 



FIG. 5. Theoretical sketch of the total magnetic field lines 
B for three diamagnetic cylinders with length-to- radius ratios 
L/R= (a) 0.2, (b) 1 and (c) 10. The lines describe the field 
direction at every point, but the field strength is not given by 
the density of lines. 



FIG. 6. Theoretical values of the surface shielding current 
K for the case L/R—l, as a function of radial and axial po- 
sitions p and z, respectively. K is normalized to the value 
Ko defined as the surface current value that would shield the 
same applied field H a for an infinite sample, that is, K~o — H a . 
The left part of the figure corresponds to the current in the 
end surfaces of the cylinder Kf aceg and the right one to the 
current in the lateral surface if lateral- The arrows mark the 
plateau of almost costant current (see text) . 



FIG. 7. Theoretical contribution from the lateral and end 
surfaces of diamagnetic cylinders to the magnetization, for 
the cases L/R=0.2, 1, and 10. Dotted, dashed and solid lines 
correspond to the contribution from the lateral surface, the 
contribution from the end surfaces, and the total, respectively. 
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